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We examine the momentum-dependent excitation spectra of indirect as well as direct resonant 
inelastic x-ray scattering (RIXS) processes in half-filled (extended) Hubbard rings. We determine the 
fundamental features of the groundstate RIXS response and discuss the experimental conditions that 
can allow for the low-energy part of these features to be distinguished in one-dimensional copper- 
oxide materials, focusing particularly on the different magnetic excitations occurring in indirect and 
direct RIXS processes. We study the dependence of spin and charge excitations on the choice of and 
detuning from resonance. Moreover, final state excitation weights are calculated as a function of 
the core-hole potential strength and lifetime. We show that these results can be used to determine 
material characteristics, such as the core-hole properties, from RIXS measurements. 
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I. INTRODUCTION 

Strongly interacting one-dimensional (ID) quantum 
systems provide the opportunity to observe distinctive 
signatures of theoretical predictions, such as the separa- 
tion of fermionic excitations into charge and spin compo- 
nents. 1 Prototypical materials which exhibit ID behavior 
are the corner-sharing copper-oxide compounds S^CuOs 
and SrCu02, in which 3D antiferromagnetic correlations 
disappear above a temperature Tjy of a few Kelvin and 
the majority of electronic properties comes from elec- 
trons moving on almost independent chains. A variety 
of experimental techniques has been used to verify the 
ID character of the excitations in these materials 2 6 and 
the existence of the expected spin-charge separation. 7 10 
These observations have induced further interest in cross- 
examining the results of different experimental probes, in 
order to test the quantitative agreement with theoretical 
descriptions of the underlying quantum mechanical pro- 
cesses. 

Resonant inelastic x-ray scattering (RIXS) is a partic- 
ularly promising technique for research in this direction, 
since it can probe multiple spin and charge excitations 
on the same footing. 11 In RIXS, incident x-ray photons 
are used that have a frequency close to the correspond- 
ing excitation energy for the promotion of a core electron 
to an empty state in or above the valence orbital of in- 
terest, the two alternatives defining direct and indirect 
RIXS respectively. 12 ' 13 In the intermediate state the re- 
sulting core hole perturbs the valence electron system lo- 
cally, until an electron decays into the empty core state, 
emitting an x-ray photon shifted in energy, momentum 
and polarization with respect to the incoming one. RIXS 
experiments at the Cu K-edge have been performed on 
the aforementioned ID cuprates 14 18 and the dominant 
spectral features have been attributed to d-d excitations. 
Spin excitations were not distinguished in the earlier ex- 
periments, because of resolution limitations of the instru- 
mentation used. These limitations can be overcome in 
today's experimental setups, as has been demonstrated 
by the observation of spin excitations in other ID 19 or 



2D 20 ' 21 cuprates and iridates, 11,22,23 although RIXS mea- 
surements for the simpler Sr-based chain cuprates are still 
awaited. 24 

The distinct form of the theoretical expression for the 
RIXS cross section provides a unique way to study prop- 
erties of models for one-dimensional condensed matter 
systems. The shortness of the intermediate state lifetime 
has been exploited to perform expansions for the mag- 
netic indirect RIXS cross section of Mott insulators. 25 27 
The case of magnetic excitations in direct RIXS has been 
investigated in this manner 28 and also by an effective 
operator method, 29 recently applied to the case of the 
quasi-lD compound TiOCl. 30 Apart from model fits to 
experimental data, 15,16 exact numerical investigations of 
indirect RIXS on small systems have been performed for 
the ID (extended) Hubbard, 31 Heisenberg, 32,33 t- J 33 and 
dp M models, whereas only a calculation of direct RIXS 
spectra for the Heisenberg and t-J models has appeared 
so far. 33 Furthermore, in the existing works, only the 
momentum dependence of RIXS spectra at specific reso- 
nance conditions has been obtained. It has been estab- 
lished experimentally that the character of excitations 
present in the final state of RIXS depends strongly on 
the detuning from resonance. Therefore, it is of interest 
to determine the effect of detuning on final state exci- 
tations in a generic model and see how it compares to 
the behavior of real systems. Moreover, the effect of the 
intermediate-state core hole, which will play an increas- 
ingly important role with the advent of time-resolved 
RIXS experiments, has been considered only heuristically 
so far. 

The purpose of this paper is to present a systematic 
numerical study of both indirect and direct RIXS pro- 
cesses in a ID Mott-insulating system. To this end, the 
Lanczos exact diagonalization method is used to calcu- 
late the RIXS spectra for the Hubbard and extended 
Hubbard models at half filling. The focus is on deter- 
mining fundamental properties and extracting trends of 
the models under changes in the core- hole parameters, 
rather than fitting experimental data. However, as the 
(extended) Hubbard model is believed to capture the es- 
sential physics of ID chain cuprates at low energies, our 
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results should be adequate to make qualitative predic- 
tions for the experimental outcome in this energy range. 
In Sec. II the theoretical formulation leading to the eval- 
uated RIXS cross section is briefly reviewed. In Sec. Ill 
the RIXS response is calculated as a function of incoming 
photon energy. Moreover, the dependence of the RIXS 
spectra upon the values of core-hole potential and life- 
time is carefully investigated and it is shown how the 
core-hole parameters can be related to direct RIXS ex- 
periments. Our conclusions are summarized in Sec. IV. 



II. RIXS CROSS SECTION 

A. Kramers - Heisenberg formula 

The unperturbed Hamiltonian Hq of choice for the TV- 
site, ID system we wish to study is the single-band Hub- 
bard model: 
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FIG. 1. (Color online) Illustration of RIXS processes in half- 
filled copper oxides, (a) Indirect RIXS process, (b) direct 
RIXS process without spin flip, (c) direct RIXS process with 
spin flip. For a spin flip to occur, the core-hole spin must also 
be flipped (dotted arrow). 



where cj a creates electrons with spin a at site j, hj a = 
t is the hopping integral and U the on-site 
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Coulomb repulsion strength. In one dimension, any non- 
zero U/t ratio is enough to guarantee a Mott-insulator 
ground state with energy E g at half filling. 1 For compar- 
ison, we also study the extended Hubbard model with an 
additional nearest-neighbor interaction 
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with V = 2t. The core hole present in the intermediate 
state acts as a localized potential scatterer for the 3d 
electrons and can be expressed as 
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where fij a is the number operator for core holes, hj a = 
1 signifying the presence of a core hole with spin a on 
site j. The core hole interacts via an attractive force 
parameterized by V c > with the total density hj = 
htf + fiji in the valence band described by Ho- The 
dipole operator that creates the intermediate state, as 
schematically illustrated in Fig. 1, is given by 
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where hj a creates a core hole and d}- a a photoexcited elec- 
tron at position Rj; k is the incoming photon momentum, 
corresponding to frequency uuk- Here we do not consider 
photon polarizations and geometrical factors. The RIXS 
cross section is given by the Kramers - Heisenberg for- 
mula: 35 
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where the sum is over all final states |/) with energies 
Ef. The RIXS amplitude has been defined as 
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where uji n = uok — ^res is the detuning of the incoming 
photon frequency from resonance and E n and Y are the 
energy and lifetime of the intermediate states. The fact 
that core holes are usually localized allows us to rewrite 
the above expression as 
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where q = k — k f is the momentum transfer and Hj = 
Ho — V c fij is a locally perturbed Hamiltonian. For indi- 
rect RIXS, one has d\ = 1, i.e., the only impact of the 
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photon is the creation of the core hole, as illustrated in 
Fig. 1(a). For direct RIXS dj a = Cj a photo-excites an 
electron into the band modelled by the (extended) Hub- 
bard model, see Fig. l(b,c). In the case of direct RIXS, 
care must be taken to sum over both the RIXS processes 
involving a photo-excited up electron and the the one in- 
volving a down electron, because the singlet character of 
the ground state is otherwise not preserved. 



B. Method 

We use the Lanczos method 36,37 to obtain the ground 
state of finite Hubbard rings and then to evaluate the 
RIXS cross section of Eq. (5). Even though the RIXS 
final state is translationally invariant, the intermediate- 
state Hamiltonian Hj, which is needed to evaluate the 
Green's functions acting on the initial state, is not. 
Translational invariance can thus not be used to reduce 
the Fock space size and all results presented here are for 
14-site rings. RIXS spectra of smaller even-site, as well 
as of 16-site rings have also been selectively calculated 
and were found to be qualitatively and quantitatively 
consistent with the 14-site results. 

Related numerical calculations of RIXS spectra using a 
similar method have been carried out 31 for a few specific 
cases and the results are in agreement with ours. We use 
t as unit of energy, the on-site Coulomb repulsion is kept 
fixed at U = 10£, and where the extended Hubbard model 
is studied, we set V = 2t. We have used Lanczos exact 
diagonalization to explore signature features of the RIXS 
response and its dependence on the remaining parameters 
of the Hamiltonian describing the core hole potential and 
the energy of the incoming photon. 



III. RESULTS 
A. Spin excitations in RIXS spectra 

The outcome of an indirect RIXS process depends 
strongly on the resonance condition. Calculation of the x- 
ray absorption (XAS) spectrum for the Hubbard model 
shows that there are two main resonances for indirect 
RIXS: 38 one corresponds to intermediate states where 
the valence band at the core-hole site is doubly occupied 
and the core hole is therefore well-screened, at energies 
~ U — 2V C with respect to the initial state, and the other 
to states where the core-hole site is singly occupied and 
the core hole is poorly screened, at energies ~ — V c with 
respect to the initial state. In states with double occu- 
pancies at the core- hole sites, exchange interactions are 
not favorable and therefore spin excitations are not pro- 
nounced, whereas they are favored when core-hole sites 
are singly occupied. This is evident in Fig. 2(a), where 
the low-energy parts of the indirect RIXS spectra taken 
close to the two resonances are presented. 
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FIG. 2. (Color online) Overview of magnetic excitations in 
calculated RIXS spectra for the Hubbard model on a half- 
filled, 14-site ring. The parameters are U/t = 10, V c /t = 15, 
uJin/t — —19, Y/t — 1, unless otherwise noted. The peaks 
are Lorentzians of width O.lt. Note that RIXS spectra in (a) 
and (c) are rescaled by the factors in parentheses. First row: 
(a) Indirect RIXS spectra for input energies corresponding to 
well-screened (solid line) and poorly screened (dashed line) 
intermediate states, (b) AS = direct RIXS spectrum (solid 
line) compared to indirect RIXS (dashed line), (c) AS = 
1 direct RIXS spectrum (solid line) compared to AS = 
direct RIXS (dashed line). Second row: Spin DSF for the 
same system (solid line) compared to AS = 1 direct RIXS 
spectrum (dashed line) for (d) Y/t = 0.2, (e) Y/t = 0.5, and 
(f) "instant" (Y/t = oo) RIXS process. 



In direct RIXS on the other hand, a double occupancy 
is created at the core-hole site by default at half filling 
and therefore there is only one resonance, corresponding 
to well-screened intermediate states. The form of the 
RIXS scattering amplitude (7) suggests that in direct 
RIXS there can be processes where a = a' (AS = 0) 
or a 7^ a' (AS = 1). In a direct RIXS experiment both 
processes take place. For a = 1 process to occur 
however, a spin-orbit interaction between the spin and 
angular degrees of freedom of the core hole is necessary 
(see Fig. 1). Hence, the weight ratio of the two direct 
RIXS processes depends heavily on the strength of the 



core-hole spin-orbit coupling. 

The low-energy parts of the spectra for the two direct 
RIXS processes are depicted in Figs. 2(b,c). It is seen 
that excitations in the low-energy part of the AS = 
process spectrum are almost identical to those occurring 
in indirect RIXS [Fig. 2(a)], although in this case they are 
caused by the decoupling of the intermediate state dou- 
blon into an antiholon and a spinon. The AS = 1 pro- 
cess spectrum is very similar to the spin dynamical struc- 
ture factor (DSF) S(q) for the studied system, shown in 
Figs. 2(d-f), and corresponds mainly to single-magnon 
excitations. Indeed, by setting r = 1/r — >> oo, i.e., by as- 
suming an instantaneous RIXS process with an infinitely 
fast decay of the core hole, S(q) is obtained exactly, see 
Fig. 2(f). The differences between S(q) and the AS = 1 
direct RIXS spectrum in Figs. 2(d,e) are thus a direct 
consequence of the core hole's finite lifetime in the inter- 
mediate state. 

The aforementioned differences between S(q) and the 
AS = 1 direct RIXS spectrum consist of changes in the 
weight distribution along the Brillouin zone, while the 
total weight remains constant, as will be shown below. 
Even though the A# = 1 direct RIXS spectrum may ap- 
pear to be almost insensitive to the value of T, the effect 
of a finite core-hole lifetime is actually important for the 
determination of the lineshape, even though the main ex- 
citation in a AS = 1 direct RIXS process - namely the 
single magnon - does not depend on the presence of the 
core hole. On the contrary, indirect and AS = 1 direct 
RIXS processes depend crucially on core-hole properties 
(see discussion below). 

The generic features of the spectra in Fig. 2 are similar 
to those found in linear spin- wave theory results, 29 but 
they have more structure. For example, in A S = 1 di- 
rect RIXS several peaks arise when moving towards the 
Brillouin zone boundary, apart from the single peak of 
single-magnon excitations. Also, the two- magnon spec- 
tral function is found to be finite at the zone boundary 
using spin- wave theory, but in Figs. 2(b,c) we see that 
the AS = direct RIXS amplitude vanishes there. This 
means that two-magnon excitations at q = ir can arise 
only in AaS = 1 direct RIXS. Our results are also con- 
sistent with numerical Bethe ansatz calculations of dy- 
namical spin and spin-exchange structure factors for the 
Heisenberg model, 32 which are the zero-lifetime limits of 
AS = 1 direct and indirect RIXS processes respectively. 

The AS = 1 part of the direct RIXS response is one or- 
der of magnitude larger than that of the AS = process. 
This is because in the AS = 1 case there will be a magnon 
excitation in the final state even if no scattering with the 
core-hole potential occurs in the intermediate state. In 
the AS = case, at least two scatterings of electrons 
off the core-hole potential are needed for magnetic exci- 
tations to arise. In direct RIXS, the incoming photon's 
angular momentum couples to the core electron's angular 
momentum, which is in turn coupled to the electron spin 
via spin-orbit coupling. This implies that it is in prin- 
ciple possible to alter the ratio of the two processes by 
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FIG. 3. Total spin excitation weight, summed over all inequiv- 
alent momenta q in the first Brillouin zone, as a function of 
cjin in (a) indirect, (b) AS = and (c) AS = 1 direct RIXS 
respectively, for the same system as in Fig. 2. The thick (thin) 
solid lines are the corresponding XAS spectra with Txas set 
to 0.5£ (O.Olt). The arrows indicate the XAS peaks at which 
the magnetic RIXS signal is maximal. 



appropriate choice of the incoming beam polarization and 
thus observe experimentally the respective change in the 
character of the spin excitations. An interesting feature 
seen in Fig. 2(c) is the difference between the one-magnon 
and the two-magnon peak locations. This difference de- 
creases with increasing chain length and becomes zero for 
long chains. The fact that the two RIXS responses scale 
differently with increasing system size is related to the 
different size of the excitations present in each case. It 
is noteworthy that experimental characterization of an- 
other quasi- ID cuprate, CaCu203, suggests the presence 
of weakly coupled, finite-size chain segments, 13-14 cop- 
per sites in length. 39 Such an energy difference might thus 
reveal signatures of chain breaks, in addition to an analy- 
sis of the static magnetic susceptibility. 40 If this is indeed 
the case, then our results show that this energy difference 
can be resolved in direct RIXS experiments and can pos- 
sibly be used to distinguish the different spin-excitation 
weights. Before proceeding, we note that all the results 
presented in Fig. 2 are qualitatively the same and quan- 
titatively very similar for the extended Hubbard model, 
for values of up to 2t for the nearest-neighbor Coulomb 
repulsion. 



B. Resonance condition and detuning 

The type of excitations present in the final state of 
RIXS is defined not only by the type of transition in the 
absorption step, but also by the position relative to the 
main resonance peak. By appropriately tuning uo- ini spe- 
cific excitations can thus be favored or hindered. This is 
evident for example in the dependence of the total spin- 
excitation weight on uo- in around the XAS peaks, as pre- 
sented in Figs. 3(a) and (b). The momentum dependence 
of the charge sector excitations can be explained within 
a free particle-hole (doublon-hole) excitation picture and 
is very similar to the one obtained in electron energy- 
loss spectroscopy (EELS) measurements. 9 The energy- 
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loss range of charge excitations in RIXS becomes nar- 
rower with increasing momentum and the main peak 
shifts to higher energy loss. 31 Similar spectral shifts of 
charge excitation peaks are observed in numerical re- 
sults for two-dimensional cuprates. 38 ' 41 A finite nearest- 
neighbor Coulomb repulsion V is equivalent to an attrac- 
tion between neighboring doublon-hole pairs and leads 
to a transfer of charge spectral weight to lower energies 
for all momenta. 31 Apart from this shift, the results pre- 
sented in Fig. 3 remain qualitatively unchanged for values 
of up to V = 2t. 

It is noted that the maxima of absorption in Figs. 3(a) 
and (b) do not coincide exactly with the maxima of spin 
excitations. Upon increasing the XAS spectrum resolu- 
tion, it is seen that the maximal RIXS spin excitation 
response lies close to specific states within the broad 
peaks. Therefore, in order to study the dependence of 
excitations on the core-hole potential strength, uji n must 
be carefully selected in each case. For this purpose, the 
cji n -dependence of the spectrum is investigated first for 
each value of V c and the value that yields the largest 
spin excitation weight is chosen as a reference. For ex- 
ample, the arrows in Fig. 3 point to the appropriate c<Ji n 
values for the case of V c = 15t. After determining the 
reference states, the ^-dependence of the total weights 
can be studied. The scaling of the total magnetic exci- 
tation weight, summed over all inequivalent momenta in 
the first Brillouin zone, is shown in Fig. 4 for indirect and 
direct RIXS processes. 



C. The role of the core-hole potential magnitude 

Figure 4 shows the relative spin excitation spectral 
weight as a function of the core-hole potential for three 
different core- hole lifetime values. As was mentioned 
above, the experimental direct RIXS spectra contain 
both the AS = and AS = 1 contributions. It is in prin- 
ciple possible to deduce the contribution of each of the 
two processes by comparing the momentum-dependence 
of the experimental spectrum lineshape (with all geomet- 
rical factors stripped off) to the lineshapes of Figs. 2(b,c). 
One can then compare with the results presented in Fig. 4 
and determine possible values for the core-hole potential 
and lifetime, which are typically considered as fit param- 
eters. We note that this estimation scheme is general 
and therefore applicable to more complicated models, 
provided that the spin excitations can be clearly distin- 
guished from other excitations, like e.g. d-d excitations 
in multiple-orbital models. 

The spin-excitation weight for indirect RIXS in 
Fig. 4(a) shows a dip around V c = U for both Hub- 
bard and extended Hubbard models. This can be in- 
tuitively understood as follows: for these values of the 
core-hole potential, the attraction exerted by the core 
hole to an electron on a different site is approximately 
canceled by the effect of the on-site Coulomb repulsion 
due to the electron already occupying the site with the 
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FIG. 4. (Color online) Total spin excitation weight percentage 
as a function of the core- hole potential in (a) indirect RIXS, 
(b) AS = direct RIXS and (b) AS = 1 direct RIXS for 
Hubbard (empty symbols) and extended Hubbard (full sym- 
bols) models, for three values of the inverse lifetime I. The 
value for the detuning oj[ n of the incoming photon energy from 
resonance at each point is defined by the procedure outlined 
in the text. All other relevant parameters are as in Fig. 2. 
(d) Core-hole lifetime dependence of magnetic excitations in 
all RIXS processes for V c = 15t. 



core hole. There is therefore no cause for scattering and 
spin exchange processes are consequently suppressed. We 
find no such dip in AS = direct RIXS, a fact which 
demonstrates the different origin of the same excitations 
in the two RIXS processes. As can be seen in Fig. 4, 
the inclusion of a nearest-neighbor Coulomb repulsion 
with V/t = 2 does not substantially alter spin-excitation 
weights. 



D. RIXS core- hole lifetime 

It has been noticed in previous studies that the value 
of the intermediate state lifetime 1/T affects the weight 
of spin excitations in model calculations. 34 The reason is 
that spin exchange processes are slow compared to charge 
excitations. This is indeed the case for indirect RIXS. 
In direct RIXS however, the single-spin flip is, as dis- 
cussed above, a zero-order process in terms of electron 
- core -hole scattering, and the total excitation weight 
does therefore not depend on cji n , r or F c , even though 
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the lineshape depends on these parameters, as is shown 
in Figs 2(d-f) for the case of T. This implies that in 
direct RIXS, evidence of single-spin excitations should 
be experimentally observable in ID chain materials, as 
is the case at the Cu L edge in 2D cuprates 20,21 and in 
Sri4Cu2404i. 19 Regarding the double spin- flip processes 
in RIXS, the core- hole lifetime has been shown to be suf- 
ficient to detect such spin excitations in the case of the 
the Cu K edge in 2D cuprates. 20 ' 42 It is thus justified to 
assume that the intermediate-state lifetime is similarly 
long enough in chain systems, suggesting that spin exci- 
tations due to exchange interactions should be present in 
experimental RIXS spectra. 

As can be seen in Fig. 4(a-c), magnetic excitation 
weights remain unchanged for large ranges of V C1 in both 
direct and indirect RIXS. Since it is very likely that 
the physical value of the core-hole attraction is within 
these ranges, the core-hole lifetime dependence of mag- 
netic excitations shown in Fig. 4(d) can be used to iden- 
tify the lifetime range corresponding to experimental ob- 
servations. The currently used estimates based on the 
XAS core- hole lifetime, as well as experimental RIXS ev- 
idence, 43 suggest that T/t < 1, a fact that underlines the 
significance of the effects of the finite core-hole lifetime 
presented above. This approach to obtain core- hole life- 
times is usable even in cases where the lifetime does not 
correspond to the natural line width. In the future, it 
may thus become particularly useful for the interpreta- 
tion of time-resolved RIXS measurements, in which the 
intermediate-state lifetime will be controlled via stimu- 
lated emission. 



IV. CONCLUSIONS 

We have determined the direct and indirect RIXS re- 
sponses of the half-filled, single-band Hubbard model in 
one dimension by using exact diagonalization, system- 
atically varying relevant parameters. In the low-energy 
sector, our results exhibit the main magnetic excitations 
allowed theoretically, namely two-magnon excitations in 



the indirect and AS = direct RIXS processes and 
single-magnon excitations in AS = 1 processes. The 
strong dependence of the total spin excitation weight 
on the resonance condition, as well as on the detuning 
from resonance, is apparent in the RIXS spectra. The 
characteristics of spin excitations in all RIXS processes 
remain qualitatively unchanged by the introduction of a 
nearest-neighbor Coulomb repulsion with magnitude of 
up to 2t. The general features of the calculated spectra 
are accessible to experiment and it is of interest to see 
to what extent the obtained generic results can describe 
measured magnetic features of ID single-chain cuprates, 
which have yet to appear. 

Furthermore, we have studied the dependence of spin 
excitations created by RIXS on the strength of the core- 
hole potential and we have presented a method for esti- 
mating the core-hole properties of interest using our or 
equivalent results, in conjunction with the respective ge- 
ometrical factors. It is also seen that spin excitations in 
indirect RIXS are suppressed for values of V c close to the 
on-site Coulomb repulsion strength U due to suppression 
of electron scattering off the core hole. The total weight 
of magnetic excitations in AS = 1 direct RIXS does not 
depend on incoming photon energy, core-hole potential 
magnitude or core- hole lifetime, and can be used as a 
reference by which the total AS = contribution can 
be distinguished in experimental data. Using the ratio 
of the magnetic weights, the core-hole properties can be 
estimated. 

During the writing of this manuscript, a preprint by 
Igarashi and Nagao appeared, 44 in which they apply their 
previously developed formalism 45 to ID systems and ob- 
tain similar results to the ones in Fig. 2. 

ACKNOWLEDGMENTS 

The authors are grateful to L. Ament, V. Bisogni and 
K. Wohlfeld for invaluable discussions. This work was 
supported by the Deutsche Forschungsgemeinschaft un- 
der the Emmy-Noether program (SK and MD). 



1 E. Lieb and F. Wu, Phys. Rev. Lett. 20, 1445 (1968). 

2 N. Motoyama, H. Eisaki, and S. Uchida, 
Phys. Rev. Lett. 76, 3212 (1996). 

3 M. Matsuda, K. Katsumata, K. M. Kojima, M. Larkin, 

G. M. Luke, J. Merrin, B. Nachumi, Y. J. Uemura, 

H. Eisaki, N. Motoyama, S. Uchida, and G. Shirane, 
Phys. Rev. B 55, R11953 (1997). 

4 A. Keren, L. P. Le, G. M. Luke, B. J. Sternlieb, 
W. D. Wu, Y. J. Uemura, S. Tajima, and S. Uchida, 
Phys. Rev. B 48, 12926 (1993). 

5 K. Kojima, M. Larkin, G. M. Luke, B. Nachumi, Y. J. Ue- 
mura, H. Eisaki, M. Motoyama, S. Uchida, B. J. Sternlieb, 
and G. Shirane, Czech. J. Phys. 46, 1945 (1996). 

6 K. M. Kojima, Y. Fudamoto, M. Larkin, G. M. Luke, 
J. Merrin, B. Nachumi, Y. J. Uemura, N. Mo- 



toyama, H. Eisaki, S. Uchida, K. Yamada, Y. En- 
doh, S. Hosoya, B. J. Sternlieb, and G. Shirane, 
Phys. Rev. Lett. 78, 1787 (1997). 

7 C. Kim, A. Y. Matsuura, Z.-X. Shen, N. Motoyama, 
H. Eisaki, S. Uchida, T. Tohyama, and S. Maekawa, 
Phys. Rev. Lett. 77, 4054 (1996). 

8 C. Kim, Z.-X. Shen, N. Motoyama, H. Eisaki, 
S. Uchida, T. Tohyama, and S. Maekawa, 
Phys. Rev. B 56, 15589 (1997). 

9 R. Neudert, M. Knupfer, M. S. Golden, J. Fink, 
W. Stephan, K. Penc, N. Motoyama, H. Eisaki, and 
S. Uchida, Phys. Rev. Lett. 81, 657 (1998). 

B. J. Kim, H. Koh, E. Rotenberg, S.-J. Oh, H. Eisaki, 
N. Motoyama, S. Uchida, T. Tohyama, S. Maekawa, Z.-X. 
Shen, and C. Kim, Nat. Phys. 2, 397 (2006). 



7 



11 L. Ament, M. van Veenendaal, T. Devereaux, J. Hill, and 
J. van den Brink, Rev. Mod. Phys. 83, 705 (2011). 

12 J. van den Brink and M. van Veenendaal, 
J. Phys. Chem. Solids 66, 2145 (2005). 

13 J. van den Brink and M. van Veenendaal, 
Europhys. Lett. 73, 121 (2006). 

14 M. Z. Hasan, P. A. Montano, E. D. Isaacs, Z.-X. Shen, 
H. Eisaki, S. K. Sinha, Z. Islam, N. Motoyama, and 
S. Uchida, Phys. Rev. Lett. 88, 177403 (2002). 

15 Y.-J. Kim, J. P. Hill, H. Benthien, F. H. L. Essler, 
E. Jeckelmann, H. S. Choi, T. W. Noh, N. Motoyama, 
K. M. Kojima, S. Uchida, D. Casa, and T. Gog, 
Phys. Rev. Lett. 92, 137402 (2004). 

16 D. Qian, Y. Li, M. Hasan, D. Casa, T. Gog, Y. Chuang, 
K. Tsutsui, T. Tohyama, S. Maekawa, and H. Eisaki, 
J. Phys. Chem. Solids 66, 2212 (2005). 

17 S. Suga, S. Imada, A. Higashiya, A. Shigemoto, S. Kasai, 
M. Sing, H. Fujiwara, A. Sekiyama, A. Yamasaki, C. Kim, 
T. Nomura, J. Igarashi, M. Yabashi, and T. Ishikawa, 
Phys. Rev. B 72, 081101 (2005). 

18 J. W. Seo, K. Yang, D. W. Lee, Y. S. Roh, J. H. Kim, 
H. Eisaki, H. Ishii, I. Jarrige, Y. Q. Cai, D. L. Feng, and 
C. Kim, Phys. Rev. B 73, 161104 (2006). 

19 J. Schlappa, T. Schmitt, F. Vernay, V. N. Strocov, 
V. Ilakovac, B. Thielemann, H. M. R0nnow, S. Van- 
ishri, A. Piazzalunga, X. Wang, L. Braicovich, G. Ghir- 
inghelli, C. Marin, J. Mesot, B. Delley, and L. Patthey, 
Phys. Rev. Lett. 103, 047401 (2009). 

20 J. P. Hill, G. Blumberg, Y.-J. Kim, D. S. Ellis, 
S. Wakimoto, R. J. Birgeneau, S. Komiya, Y. Ando, 
B. Liang, R. L. Greene, D. Casa, and T. Gog, 
Phys. Rev. Lett. 100, 097001 (2008). 

21 L. Braicovich, J. van den Brink, V. Bisogni, M. M. Sala, 
L. J. P. Ament, N. B. Brookes, G. M. De Luca, M. Sal- 
luzzo, T. Schmitt, V. N. Strocov, and G. Ghiringhelli, 
Phys. Rev. Lett. 104, 077002 (2010). 

22 J. Kim, D. Casa, M. H. Upton, T. Gog, Y.-J. 
Kim, J. F. Mitchell, M. V. Veenendaal, M. Daghofer, 
J. van den Brink, G. Khaliullin, and B. J. Kim, 
arXiv: 11 10:0759 (2011). 

23 M. van Veenendaal, arXiv: 1106.0640 (2011). 

24 J. Schlappa, K. Wohlfeld, K. J. Zhou, M. Mourigal, 
M. W. Haverkort, V. N. Strocov, L. Hozoi, C. Mon- 
ney, S. Nishimoto, S. Singh, A. Revcolevschi, J. S. Caux, 



L. Patthey, H. M. R0nnow, J. van den Brink, and 
T. Schmitt, submitted to Nature (2011). 

25 J. van den Brink, Europhys. Lett. 80, 47003 (2007). 

26 L. J. P. Ament, F. Forte, and J. van den Brink, 
Phys. Rev. B 75, 115118 (2007). 

27 T. Nagao and J. I. Igarashi, 
Phys. Rev. B 75, 214414 (2007). 

28 L. J. P. Ament, G. Ghiringhelli, M. M. 
Sala, L. Braicovich, and J. van den Brink, 
Phys. Rev. Lett. 103, 117003 (2009). 

29 M. W. Haverkort, Phys. Rev. Lett. 105, 167404 (2010). 

30 S. Glawion, J. Heidler, M. W. Haverkort, L. C. Duda, T. 
Schmitt, V. N. Strocov, C. Monney, K. J. Zhou, A. Ruff, 
M. Sing, and R. Claessen, Phys. Rev. Lett. 107, 107402 
(2011). 

31 K. Tsutsui, T. Tohyama, and S. Maekawa, 
Phys. Rev. B 61, 7180 (2000). 

32 A. Klauser, J. Mossel, J.-S. Caux, and J. van den Brink, 
Phys. Rev. Lett. 106, 157205 (2011). 

33 F. Forte, M. Cuoco, C. Noce, and J. van den Brink, 
Phys. Rev. B 83, 245133 (2011). 

34 K. Okada and A. Kotani, 



J. Phys. Soc. Jpn. 75, 044702 (2006). 



H. 



A. Kramers and W. Heisenberg, 



Z. Phys. 31, 681 (1925). 

36 C. Lanczos, J. Res. Natl. Bur. Stand. 45, 255 (1950). 

37 C. Lanczos, J. Res. Natl. Bur. Stand. 49, 33 (1952). 

38 K. Tsutsui, T. Tohyama, and S. Maekawa, 
Phys. Rev. Lett. 83, 3705 (1999). 

39 V. Kiryukhin, Y. J. Kim, K. J. Thomas, F. C. Chou, 
R. W. Erwin, Q. Huang, M. A. Kastner, and R. J. Birge- 
neau, Phys. Rev. B 63, 144418 (2001). 

40 J. Sirker, N. Laflorencie, S. Fujimoto, S. Eggert, and I. Af- 
fleck, Phys. Rev. Lett. 98, 137205 (2007). 

41 C. J. Jia, C.-C. Chen, A. P. Sorini, B. Moritz, and 
T. P. Deveraux, arXiv: 10 10.3446 (2011). 

42 V. Bisogni, L. Simonelli, L. J. P. Ament, F. Forte, 
M. M. Sala, M. Minola, S. Huotari, J. van den Brink, 
G. Ghiringhelli, N. B. Brookes, and L. Braicovich, 
arXiv: 1010.4725 (2010). 

43 V. Bisogni, private communication. 

44 J. Igarashi and T. Nagao, arXiv:l 11 1.5078 (2011). 

45 J. Igarashi and T. Nagao, arXiv: 1104.4683 (2011). 



